Localization and derealization of two-dimensional discrete solitons pinned to linear 
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We study the dynamics of two-dimensional (2D) localized modes in the nonlinear lattice described 
by the discrete nonlinear Schrodinger (DNLS) equation, including a local linear or nonlinear defect. 
Discrete solitons pinned to the defects are investigated by means of the numerical continuation from 
the anti-continuum limit and also using the variational approximation (VA), which features a good 
agreement for strongly localized modes. The models with the time-modulated strengths of the linear 
or nonlinear defect are considered too. In that case, one can temporarily shift the critical norm, 
below which localized 2D modes cannot exists, to a level above the norm of the given soliton, which 
triggers the irreversible derealization transition. 

PACS numbers: 03.75.Lm, 03.75.Kk, 03.75.-b 



I. INTRODUCTION 

The discrete nonlinear Schrodinger (DNLS) equations 
constitute a universal class of models which are pro- 
foundly interesting in their own right, as dynamical sys- 
tems, and serve as models of a number of physical sys- 
tems in nonlinear optics 1| , studies of Bose-Einstein con- 
densates 0,0], and in other contexts In particular, 
soliton solutions to the DNLS equation in one, two, and 
three dimensions represent fundamental dynamically lo- 
calized modes in discrete media. Experimentally, one- 
and two-dimensional (ID and 2D) discrete solitons have 
been created in nonlinear optical systems of several types 

a- 

An important ingredient of DNLS models is repre- 
sented by local defects. They are interesting as addi- 
tional elements of the lattices [4], and find important 
physical realizations. In particular, they may describe 
various strongly localized structures in photonic crystals 

5 1 , including nanocavities @ , as well as micro-resonators 

7|, and quantum dots @. 

The objective of the present work is to consider the 
interaction of 2D discrete solitons with local linear and 
nonlinear defects in DNLS lattices. The defect may be 
concentrated at a single site of the lattice, or it may be 
shaped as a Gaussian of a finite width. After introducing 
the model in Section 2, we consider stationary solitons 
pinned by the defects in Section 3. The analysis is per- 
formed by means of the variational approximation (VA) 
and numerical methods. In particular, the pinned-soliton 
families feature a specific combination of stable and un- 
stable parts. In Section 4, we consider possibilities of 



the application of the "management" [9( to 2D discrete 
solitons, using the linear and nonlinear defects whose am- 
plitudes slowly vary in time. In this direction, we ana- 
lyze a possibility to trigger a derealization transition by 
means of this method, which may be used in the design of 
switching schemes in photonics. Previously, the induced 
transition to the derealization was demonstrates in uni- 
form lattices with the inter-site coupling strength subject 
to the time modulation The paper is concluded by 
Section 5. 



II. THE MODEL 

We consider the following model based on the 2D 
DNLS equation with a local defect: 



+ (a + W 1l:m )\ 



o. 



(i) 



where the overdot stands for the time derivative, 

^2^n,m = ^?i,m+l "h ^n.m— 1 "h ~T" Hn—l.m ^^n 7 m 

is the 2D discrete Laplacian, the coupling constant of 
the lattice will be fixed by scaling, J = 1, unless an- 
other choice of J is specified explicitly, and a = +1 and 
— 1 corresponds to the attractive and repulsive nonlin- 
earities, respectively (a — for the linear lattice). Fur- 
ther, functions V n _ m and W n ^ m , which account for the 
linear and nonlinear components of the defect, are taken 
as Gaussians profiles, 



Ki,m = v exp [-{n 2 + m 2 )/A v ] , 
W n , m = w exp [~(n 2 + m 2 )/A w ] , 



(2) 
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with respective strengths v, w and widths A„, A„. In this 
notation, positive and negative strengths correspond to 
the attractive and repulsive defect, respectively. In fact, 



2 



we will consider the linear and nonlinear defects sepa- 
rately. Models of 2D nonlinear lattices with other types 
of local defects were considered earlierjljj , including de- 
fects induced by edges of the lattice jl2l ]. 
Looking for the stationary solutions, 



= U n , m exp(-iwt), 



we arrive at the nonlinear eigenvalue problem, 



(3) 



Equation {5} can be derived from the Lagrangian, 

L = ^ ^ t ( U m,n'i 1 m,n ~ u m,nU mn ) 



+ W m) „ (u m+l n + U m n+1 ^j 4:\u m . n \ 



2^ 



u mo,o 



(6) 



+ (T + VK n , m )|C/„ im | 2 C/„ 



= 0, 



(4) 



for real frequency w and the profile of the stationary 
discrete mode U m<rl , which may be complex, in the 
general case. In the absence of the defect, Eq. Q 
gives rise to the linear dispersion relation featuring the 
phonon band, above and below which one can find non- 
linear modes, described by respective curves N(lu), with 
TV = J2 m n \ u m,n\ 2 being the norm (power) of the nonlin- 
ear state. The important difference between 2D and ID 
settings is that, in the latter case, the fundamental single- 
peak mode (the discrete soliton of the Sievers-Takeno 
type) persists in the limit of N — > 0, while all the 2D 
solitons are bounded by a critical norm, N cr , below which 
localized modes do not exist Accordingly, an abrupt 
derealization (decay) of discrete 2D solitons was pre- 
dicted in the case when the inter-site coupling constant 
would exceed a certain critical value [Toj |. In the follow- 
ing we show that, introducing the defect with the linear 
and nonlinear components and varying their strengths, 
or the lattice coupling constant, one can change the crit- 
ical norm, iV cr , and thus control the transition to the 
derealization. 



III. STATIONARY DISCRETE SOLITONS 
PINNED TO THE DEFECT: THE VARIATIONAL 
APPROACH AND NUMERICAL RESULTS 

A. The variational approximation 

The variational approach (VA) was successfully used 
for the study of 2D discrete solitons in uniform (defect- 
free) DNLS lattices [TJ, [Hj]. Here, we start with the 
application of the VA to the 2D lattice in the presence of 
the (5-defect localized at the origin: 



+ (v + w\u m , n \ 2 )5 m . 8 n fiU min - 



0. 



(5) 



cf. the more general form of the defect in Eq. (JSJ). Gen- 
erally, the strength of the defect may be time-dependent, 
v = v(t), w = w(t). 



To apply the VA, we adopt the following ansatz (previ- 
ously, it was used for the analysis of discrete solitons in 
the 2D DNLS equation with the cubic-quintic nonlinear- 
ity Q3): 



B. 



if m = n = 0: 



^ e -o(|m|+|n|) ) otherwise. 



(7) 



Substituting ansatz ([7]) into Lagrangian (J5J), one can 
perform the respective calculations and find the effec- 
tive Lagrangian as a sum of four terms which depend 
on three dynamical parameters, (A, B, a) (which may be 
functions of time) and represent, respectively, the kinetic 
part, on-site self-interaction, inter-site couplings, and ef- 
fects induced by the defect: 



Lkin — 



L — Lkin + L coup i + Li n t + L se if \ (8) 

2ie~ 2a 



(1 - e- 2a f 



A* A- AA* 



- ( B*B - BB" 



Lseif = 2a\A\ 4 



(1 - e- Aa ) 



Lcoupi — 8|A|' 



(1 - e- 2a Y 



(3e~ a - 



-4(AB* +A*B)e- a - A\B\ 
Ldofcct = v\B\ 2 + ^IBI 4 . 



(9) 
(10) 

(11) 
(12) 
(13) 



The substitution of the harmonic time dependence, 
A(t) = Ae~ %ut , B{t) = Ber^\ where A and B are real 
amplitudes, casts the Lagrangian into the following sta- 
tionary form, with C = e~ a : 



A 2 C 2 



Lkin = 4w 
L S cif = 2a A 4 
-^coupi = 8A 2 



(i-c 2 y 

c 4 
(i-c 4 ) 
c 2 



ujB 2 



2 2 



(14) 
(15) 



(i-c 2 y 



(3C - C* 3 - 2) (16) 



+8 ABC - 4B 2 , 



L, 



defect 



= vB z 



-B 4 



(17) 
(18) 



Fixing frequency uj and numerically solving the ensu- 
ing system of the Euler-Lagrange equations, which fol- 
lows from the effective Lagrangian, dL/dA = dL/dB = 



3 



dL/dC = 0, one can find the set of variational parame- 
ters ( A, B, a) and the norm corresponding to ansatz ([7]), 

^ = S 2 ,5 m ,o^,o + A 2 ^(l-W Il ,o)C 2 (l m l+l"l) . (19) 



B. The linear defect 

Curves N(u>) for the discrete solitons produced by the 
VA, along with their counterparts, obtained by dint of the 
direct numerical solution based on the continuation from 
the anti-continuum limit (see Appendix A), are displayed 
in Fig. [TJ for the pure linear defect. 
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FIG. 1. (Color online) (a) Comparison of the dependence 
of the norm on the frequency for the fundamental discrete 
solitons, as predicted by the VA, with the use of Eq. (|19|) 
(dotted blue lines), and its counterpart produced by the nu- 
merical solution of the stationary version of Eq. ([5]) (solid 
black and dashed red lines), for three different strengths of 
the linear defect (w = 0): v = —1; 0; +1 (from top to bottom; 
recall that v > and v < correspond to the attractive and 
repulsive defects, respectively). The numerical results were 
obtained for a narrow defect, with width A„ = 0.1 in Eq. 
The nonlinearity is attractive (a = +1). (b) Zoom of the plot 
in (a), corresponding to the VA prediction for the norm given 
by Eq. (|19[1 . in the vicinity of u = 0. In this and all other 
figures, solid and dashed portions of the numerically gener- 
ated curves depict subfamilies of stable and unstable solitons, 
respectively. 

Figure [TJa) shows that, in the case of A„ = 0.1, which 
is practically tantamount to the defect localized at the 
single site, cf. Eq. ((5|), the VA-generated norm is in a 
good agreement with with the numerical results far from 
the limit of oj — » 0. At small to, the soliton spreads out, 
and its actual shape deviates from the exponential ansatz 
([7|) , which leads to a discrepancy between the variational 
and numerical results, as seen in Fig. [Ha). Also due 
to restrictions on the size of the domain of calculation 
(in our case we used 41 x 41 grid points) the numerical 
curves cannot be continued to extremely small values of 
oj, therefore they are not presented in Fig. Hfb). 

Naturally, the increase of the width of the defect to 
A„ = 1, which makes the shape of defect ([2]) essentially 
different from that in Eq. ([5]), leads to a stronger dis- 
crepancy in the limit oj — > 0, as seen in Fig. [5J 
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FIG. 2. (Color online) The same as in Fig. [TJa), but with 
the norm found numerically for the broad defect with width 
A„ = 1. 



We have checked the linear stability of the solutions 
along their existence curves (see details in Appendix B). 
The results are shown by means of solid and dashed por- 
tions in the figures displaying the respective numerically 
generated N(oj) curves. Note that, here and in similar 
plots displayed below, the soliton families may contain 
two distinct stability segments separated by an instabil- 
ity region. It should be noted found that, far from the 
limit of oj = 0, the Vakhitov-Kolokolov (VK) criterion, 
dN/duj < [H , correctly predicts the transition from the 
stability to instability at points where dN/doj changes 
its sign from negative to positive. However, at small oj, 
the formally applied VK criterion only partially complies 
with the linear-stability results, which may be an effect of 
boundary conditions on properties of very broad modes 
corresponding to small oj. 

It is worthy to note that, as one can conclude from 
the comparison of Figs. QJa) and El the increase of the 
width of the attractive linear defect leads to a significant 
reduction of the instability region. This trend is also 
corroborated by the VK criterion. 



C. The nonlinear defect 

Similar results were obtained for the pure nonlinear 
defect, i.e., with v = in Eqs. ([2]) and ([5]), as shown in 
Figs. [3] and [4] In particular, analyzing the results for 
the linear and nonlinear defects with different widths, we 
have found that there is a particular intermediate value 
of the width in the interval of A„.„ C (0.1, 1), at which 
the discrepancy between the numerical and variational 
solutions attains a minimum (in particular, it is A„, « 0.7 
for the nonlinear defect, see Fig [5]). 

For the sake of comparison, the VA-predicted and nu- 
merical curves N(oj) are also shown in Fig. [6] for the 
discrete solitons in the defect-free 2D lattice (of course, 
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FIG. 3. (Color online) The same as in Fig. [T] but for the pure 
nonlinear defect (v — 0), at different values of its strength: 
w = -0.5; -0.2; 0; +0.5 (from top to bottom). The width of 
the defect corresponding to the numerical curves is A w — 0.1. 
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FIG. 4. (Color online) The same as in Fig. UK a), but for the 
defect's width A w — 1. 



these results are not different from those reported in ear- 
lier works which were dealing with the 2D DNLS equation 
without defects [1, HH). Comparing these to Figs. [T]and 
121 we conclude that the discrepancy between the VA and 
numerical findings is actually smaller in the presence of 
the linear or nonlinear defect. 



Finally, it is also interesting to consider the case of the 
linear lattice (a ~ 0), with all the nonlinearity concen- 
trated only in the form of a narrow defect © with v = 
and w > 0. In Fig. [7] , the corresponding dependence 
N(lu) is shown for A w = 0.1, which implies that the 
nonlinearity is actually concentrated at the single site, 
n = m = 0. 




FIG. 5. (Color online) The comparison between the vari- 
ational and numerical results for the pure nonlinear defect 
with amplitude w = —0.5 and width A m — 0.7. In this case, 
the overall discrepancy between the VA-predicted and numer- 
ically found curves [which are defined as in Fig. [3K a)] attains 
its minimum. 




FIG. 6. (Color online) The same as in Figs. OJa) and^a), 
but for the fundamental discrete solitons in the 2D lattice 
without defects (v = w = 0), at three different values of the 
coupling constant in Eq. JT}: J — 1.2; 1; 0.5 (from top to 
bottom). 



IV. LOCALIZATION-DELOCALIZATION 
TRANSITION CONTROLLED BY TEMPORAL 
MODULATIONS OF THE LINEAR AND 
NONLINEAR DEFECT 

A. The variation of the linear defect 

In this section, we consider effects of the hysteresis for 
the discrete solitons, induced by the adiabatic variation 
of the strength of the linear defect, v(t), which eventually 
returns to the initial value, Vi = v(t = 0), cf. Ref. (l5j |: 



v(t)=vf + (vi-v f )\l-2t/tf\. 



(20) 
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FIG. 7. (Color online) The comparison of the numerically 
found dependence N(u>) in the linear lattice (a — 0), with 
the nonlinearity represented solely by the narrow nonlinear 
defect (v = 0, A w ~ 0.1)-red lines-with its counterpart in the 
usual defect-free DNLS lattice (black solid and dashed lines). 



Here Vf = v(t — tf/2) corresponds to the turning point, 
t = tt/2, while the return time is t = it. Note that the 
total norm remains a conserved quantity in the DNLS 
equation with the time-dependent strength of the local 
defect. 

First, we consider the soliton for Vi = 0, with the norm, 
N s (point A in Fig. |8]), which precisely coincides with the 
minimum of the existence curve, but for v = — 1 (point 
B). This means that, by taking Vi — and Vf = —1 in Eq. 
(|2T))) at t — tf/2, we may transfer the soliton from point 
A to B. As long as the soliton stays on the existence curve 
between these two points, it remains localized, featuring 
some changes of the profile, see details of dynamics in 
Fig. |9] (this is also valid for the dynamics between points 
A and C). However, if, in the framework of the same sce- 
nario, the amplitude of the defect falls below the critical 
value, Vf < v CI , which in the present case is t; cr = — 1, 
then the minimum of the existence curve turns out to 
be higher than the norm of the evolving mode, which is 
expected to trigger the delocalization. 

These assumptions were verified by means of the direct 
integration of Eq. ([T]) with the strength of the linear 
defect v varying in time according to Eq. (f2U)) . The 
result is presented in Fig. |H1 where the evolution of the 
density at the origin, |uo,o(i)| 2 > is shown for three values 
of the Vf, along with the corresponding initial and final 
profiles. 



B. The variation of the nonlinear defect 

Here we aim to consider effects produced by the adi- 
abatic variation of the strength of the nonlinear defect, 
following the same scenario as in Eq. (|2T)1) : 
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FIG. 8. (Color online) The illustration of the transforma- 
tion of the discrete soliton under the action of the adiabatic 
variation of the strength of the linear defect. 
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w{t)=wf + (w i -w f )\l-2t/tf\, 



(21) 
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FIG. 9. (Color online) (a) The evolution of the density at 
the origin, |«o,o(t)| 2 , corresponding to the time-modulation 
scenarios displayed in Fig. [8] The dark solid and blue dashed- 
dotted lines stand for the scenarios A-B-A (with Vf = — 1) and 
A-C-A (with Vf = 1), which are shown in panel (b) and in 
Fig. [S] The red dashed line corresponds to the case when the 
strength v(t) temporarily falls to a value below the critical 
value, Vf = —1.1 < Vcr, which leads to the delocalization. In 
the right panels, initial (top) and final density profiles with 
Vf = —1 (middle) and Vf = —1.1 (bottom) are displayed. 



with Wi — w(t = 0) and Wf = w(t = tf/2) corresponding 
to the initial strength and its value at the turning point, 
t = tf/2, with the return to the initial value at t = tj. 
As in the previous case, the norm of the solution remains 
conserved in the course of the evolution. 

Again we start with the solution at w% = (point A in 
Fig. UOp with the norm which coincides with the mini- 
mum of the existence curve for w = —0.2 (point B). This 
means that, by taking Wi — and Wf = —0.2 in Eq. (|2"Tj) 
at time t — tf/2, we will transfer the soliton from point 
A to point B, see Fig. [TUJ If the soliton stays on the exis- 
tence curve between these two points, the solution always 
remains localized (see the dynamical picture in Fig. 1111) . 
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However, if, for the same initial condition, the strength of 
the nonlinear defect is allowed to drop below the critical 
value, Wf < w CI , which in the present case is w CI = —0.2, 
then the minimum of the existence curve becomes higher 
than the norm of the solution, which again is expected 
to trigger the transition to the derealization. 
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FIG. 10. (Color online) The illustration of the soliton trans- 
formation under the action of the adiabatic variation of the 
strength of the nonlinear defect. 

These assumption were verified through the direct in- 
tegration of Eq. ([1]) with the strength of the nonlinear 
defect, w, varying in time according to Eq. (|2ip . The re- 
sults are presented in Fig. QTJ where the evolution of the 
density at the origin, \u _ (t)\ 2 , is shown for three values 

Of Wf. 
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Schrodinger) equation, which includes the local linear or 
nonlinear defect. The solitons trapped around the de- 
fects were investigated using both the VA (variational 
approximation) and numerical methods, the VA show- 
ing good agreement with the numerical findings for suffi- 
ciently narrow solitons. Then, the analysis was extended 
to the model with the linear and nonlinear defects whose 
strength was subject to the slow variation in time. In the 
latter case, one of the possibilities is the controlled onset 
of the transition to derealization. 

The analysis reported in this work can be extended by 
considering combinations of linear and nonlinear defects, 
searching for symmetric, antisymmetric, and asymmet- 
ric modes trapped by symmetric pairs of defects, and 
analyzing other "management" schemes, with the defect 
strength subject to a time-periodic modulation. 
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APPENDIX A 

Here we discuss some details on the continuation of so- 
lutions from the anti-continuum (AC) limit, which cor- 
responds to the uncoupled lattice described by Eq. ([1]) 
and (J5J) with J = [4]. Looking for stationary modes in 
the form of Eq. ([3]), in the AC limit one obtains obvi- 
ous exact solution for the one-site fundamental discrete 
soliton, 



Un 



±A, if m = n = 0; 
0, otherwise. 



(22) 



where A = \J — (uj + v) / (a + w), see Ref. [19j and ref- 
erences therein. Then, by using the standard Newton- 
Raphson method, we gradually increase the coupling con- 
stant from J = to J = 1, also increasing the width of 
the defect if needed [to proceed from the narrow defect 
§5§ to the broader one ©]. After that, by changing fre- 
quency uj of the solution, we obtain the dependence of 
norm TV on uj, and also analyze the linear stability of the 
corresponding solutions, see the next Appendix. 



FIG. 11. (Color online) The same as in Fig. [9] but for the 
temporal modulation of the nonlinear defect. 



APPENDIX B 



V. CONCLUSION 

In this work, we have considered the static and dy- 
namical properties of 2D discrete solitons in the non- 
linear lattice described by the DNLS (discrete nonlinear 



To analyze the stability of solitons, we consider per- 
turbed solutions, 

Um,„(t) = [U m , n + a m , n (t)] e^' , (23) 

Representing the small perturbations as a m>n (t) — 
{ctm,n + ifim,n) e At , substituting this into Eq. ([I} and 
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separating real and imaginary parts, we derive the fol- 
lowing system, 




where A is the stability eigenvalue, and operators L\ and 
L3 are 

L x = A 2 + V m , n + + W m , n ) {U m , n f , (25) 
L 3 = A 2 + V m , n -uj + 3(a + W m>n ) (U m>n f . (26) 



Multiplying the first equation in system (|24|) by —A and 
than using the second equation, one arrives at the eigen- 
value problem, 

Aa = L±L 3 a, (27) 

where A = —A 2 . The underlying stationary solution, 
U mn , is unstable if the spectrum of A includes complex 
or negative real values. 
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